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ABSTRACT 

We have performed a cosmological numerical simulation of primordial baryonic gas 
collapsing onto a 3 x lO^M© dark matter (DM) halo. We show that the large scale 
baryonic accretion process and the merger of few ~ 10^ DM halos, triggered by the 
gravitational potential of the biggest halo, is enough to create super sonic {M. > 10) 
shocks and develop a turbulent environment. In this scenario the post shocked regions 
are able to produced both H2 and HD molecules very efficiently reaching maximum 
abundances of nn^ ~ 10~^nH and uhd ^ few x 10~^nu, enough to cool the gas below 
lOOK in some regions. The kinetic energy spectrum of the turbulent primordial gas 
is close to a Burgers spectrum, Ek oc k~'^, which could favor the formation of low 
mass primordial stars. The solenoidal to total kinetic energy ratio is 0.65 < < 0.7 
for a wide range of wave numbers; this value is close to Rk « 2/3 natural equipartition 
energy value of a random turbulent flow. In this way turbulence and molecular cooling 
seem to work together in order to produce potential star formation regions of cold and 
dense gas in primordial environments. We conclude that both the mergers and the 
collapse process onto the main DM halo provide enough energy to develop super sonic 
turbulence which favor the molecular coolants formation: this mechanism, which could 
be universal and the main route toward formation of the first galaxies, is able to create 
potential star forming regions at high redshift. 

Key words: galaxies: formation — large-scale structure of the universe — stars: 
formation — turbulence. 



1 INTRODUCTION 

In the currently accepted paradigm to describe the universe, 
LCDM, the formation of bound dark matter structures takes 
place via hierarchical growth. While this build up is driven 
by dark matter accretion onto an already virialized dark 
matter halo, the process of merging is equally important. 
The effect of dark matter mergers in building up the dark 
matter halo has been studied in enormous detail both nu- 
merically and analytically. Due to new hydrodynamical sim- 
ulations that have achieved enough resolution, the effect of 
merging of baryons is being taken into account, including 
also the important process of transfer of baryonic matter 
between dark matter merging halos. 

One aspect of merging that has not been studied in 
detail is the development of (supersonic) turbulence in the 
merging halos due entirely to the dark matter merging pro- 
cess and its consequent effect on the environment. The 
current paradigm to power turbulence in the interstellar 
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medium is mainly through the process of SN explosions 
( Norman fc Ferrara||1996 1. However, it is possible that the 
dark matter merging process itself does power the initial 
phase of turbulence and thus star formation. 

Using numerical simulations from cosmological initial 



condition Wise & Abel ( 2007 1 and Greif et al. ( 2008 1 have 



discussed the generation of turbulent motions in primor- 
dial gas through the virialization process of halos of mass 
^ 10*Mq and ^ 5 x lO'^M©, respectively. They have shown 
that supersonic turbulent motions, which partially ionize the 
primordial gas allow an efficient formation of H2 and HD 
molecules and efficient cooling to a gas temperature below 
~ lOOK in some regions. They have argued that these turbu- 
lent low temperature regions could be sites of efficient star 
formation. 



The result of Greif et al. ( 2008 1 is somehow different 
from the conclusions of |Johnson fc Bromrn] |2006) who have 
analytically shown that in order to ionize the primordial 
gas through the DM virialization process a halo mass above 
lO^M© is needed. This apparent contradiction is easily un- 



derstood because Johnson & Bromm ( 2006 1 took into ac- 
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count the fixed halo average circular velocity for the veloc- 
ity shock waves instead a velocity distribution for primordial 
gas (which can take values higher and lower than the aver- 
age circular velocity). Summarizing, Greif et al. (20081 have 
shown that the DM virialization process of halos with mass 
> IO^Mq can generate a supersonic turbulent media (trig- 
gered by both hot and cold gas accretion) where coolants 
(H2 and HD molecules) are generated efficiently. However, 
despite of the interesting insights on the development of tur- 
bulence in cosmological simulations of primordial gas made 
by the simulation mentioned above, their authors have not 
pay attention on the impact of the halo merger process as 
a turbulence inductor and they have related the onset of 
turbulence just with the gas accretion onto the DM halos. 



On the other hand, Wise & Abel ( 2007 1 mentioned that 



the minor and major mergers are able to produce turbulence 
by Kelvin-Helmholtz instabilities, as well as by gas accre- 
tion. Despite of this mention, the both studies mentioned 
above neither characterize the turbulence with power spec- 
trum in protogalaxy environments nor characterize the gas 
conditions with probability distribution functions in order 
to have a clear idea about the primordial gas conditions in- 
side the first protogalaxies. No previous study has ever done 
this in such kind of environments, therefore we believe our 
study brings many new results to the field. 

In the last years there have been a number of works 
studying the development of turbulence under isothermal 
and non isothermal conditions. A very interesting aspect of 
supersonic gas turbulence is the effect of the gas dynamic on 
the gas chemistry (mainly on the formation of gas coolants). 
[Milosavljevic et al.| ( |201l"| have shown that compressive forc- 
ing turbulence produces H2 formation faster than in the 
solenoidal forcing case. This behavior is explained by these 
authors due to the broader compressive forcing density PDF 
which creates more £[2 molecules at initial times. The results 



of Milosavljevic et al. (2011 \ could be placed in a cosmolog- 



ical context and argue that if a DM halo is able to produce 
enough compressive supersonic turbulent motions it is pos- 
sible to produce a large amount of molecules and cool down 
the primordial gas very efficiently. 

In this paper we investigate in detail how turbulence 
is generated by both mini DM halo merging and gas col- 
lapse using a cosmological size simulation. We characterize 
the gas physical conditions computing probability distribu- 
tion functions of gas (density, temperature and velocity) and 
the gas velocity power spectrum of this cosmological simu- 
lation, which is an study not done and can be useful for 
idealized primordial star formation simulations. We follow 
the merging process in detail and investigate the onset of a 
turbulent medium inside the merging halos, including very 
small mass halos (Af ~ 10^ Mq). The most interesting find- 
ings of our work are that merging does indeed power turbu- 
lence, but most interesting is that it enhances significantly 
the coolants in the primordial gas, thus significantly lower- 
ing the temperature of the primordial gas, even below 100 
K in some regions. Furthermore, we found that the velocity 
power spectrum is nearby a Burgers spectrum oc and 
the solenoidal to total kinetic energy ratio is nearby 2/3 as 
found in idealized random fluids simulations. 

The merger inducing turbulence process can be univer- 
sal and can also take place at lower redshift, thus powering 
the initial phases of star formation in merging halos. The pa- 



per is structured as follows: in § 2 we describe the methodol- 
ogy and follow in § 3 with our main results and analysis. We 
continue with a discussion of our results in § 4 and conclude 
and summarize in 5 5. 



2 METHODOLOGY 

We used the adaptive mesh refinement (AMR) code RAM- 
SES ( |Teyssier||2002| with a modified non-equilibrium cool- 
ing module with 21 species (including H2, HD, and LiH 
molecules with their cooling functions) in order to follow 
the chemo-thermal evolution of primordial gas and to study 
the gas gravitational collapse process in pristine environ- 



ments (Prieto, Infante, & Jimenez 20081. All chemical re- 



action rates were taken from Standi et al. (19961, Galli & 



Palla (19981, and Glover & Abel (20081 



We worked in a concordance ACDM cosmological 
model: h = 0.72, f^A = 0.73, D.^ = 0.27, fib = 0.04, erg = 0.9 
and Us — 0.95. The dynamical initial conditions were taken 
from mpgrafic ( jPrunet et al.||2008[) and th e initial chemical 
abundances were taken from Galli & Palla ( 1998 1 at z = 120 



(the initial redshift of the simulation). 

We first simulated a volume of (1 Mpc/h) ^ of only DM, 
with 256^ particles. Using the HOP algorithm (Eisenstein & 
Hut|l998| ), we identified a 3 x lO'^ M© halo at z = 10. In or- 
der to better resolve the formation of that halo and to study 
its baryon component, we re-simulated the same (IMpc/h)'^ 
volume with 512'^ particles (with a particle mass of approx- 
imately ~800Mq), including gas with non-equilibrium pri- 
mordial chemistry from the beginning of the simulation at 
z=120. The gas dynamics was computed on a root grid of 
64^ computational zones, which was geometrically refined 
towards the center, increasing the spatial resolution by a 
factor of two within the central 1/8 of the volume. This 
refinement criterion was applied three times, generating 4 
nested meshes (including the root grid), each with 64^ com- 
putational cells. 

Inside the innermost 64^ mesh, the gas dynamics was 
computed using seven extra levels of adaptive mesh refine- 
ment by a factor of two. This adaptive refinement was based 
on four different criteria: i) Lagrangian refinement based on 
the number density of DM particles (a finer refinement level 
is created in cells containing more than 4 DM particles), 
ii) Lagrangian refinement based on the baryonic mass den- 
sity, iii) refinement based on the gas pressure gradient (for 
Ap/p > 2), and iv) refinement based on the Jeans' length, 
to satisfy Truelove's condition (Truelove et al. 19971. The 



pressure gradient criterion was included in order to better re- 



solve the turbulent ffow, as discussed in Kritsuk et al. ( 2006 1 



The three geometrical refinement levels, plus the additional 
seven adaptive refinement levels, give an effective spatial res- 
olution corresponding to that achieved by a uniform mesh 
of 65, 536'^ computational elements, and corresponding to a 
proper size of «1.95pc a.t z ~ 10.0. 

This setup was used from the beginning of the simu- 
lation, at 2 = 120, until z = 50. At z = 50, another geo- 
metrical refinement criterion was added (while maintaining 
the others): Uniform spatial resolution till level 15 was im- 
posed around the densest central region of the halo, creating 
a uniform mesh of 512^ elements which was refined with a 
J = 0.125 Jeans criteria till level 16 (due to technical rea- 
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sons the uniform refinement criteria was imposed gradually 
and it was completed at z=30). This high refinement region, 
covering a proper size of ~2.0 kpc at z « 10.0, served the 
purpose of better resolving the generation of turbulent mo- 
tions in the central region of the halo. The results presented 
in this work are based on the analysis of this central 2.0kpc 
region at the end of the simulation a.t z ~ 10.0. 



3 RESULTS AND ANALYSIS 

The simulation presented in this work follows the baryonic 
matter accretion process onto a « 3 x lO^M© (hereafter the 
main halo) with Tvir « 7770 K and Kirc ~ 14.6 km/s. 

At z « 14 the main halo is crossed by two interacting 
DM minihalos of ~ IO^Mq. These two interacting miniha- 
los merge and create a spinning baryonic over-density going 
out the main halo central region. Due to this violent merging 
process the gas develops a shock wave which increases the 
free electron fraction. The post shocked gas develops tur- 
bulence and creates H2 and HD molecules very efficiently 
behind the merged over-density. 

After the merging process the main halo continues ac- 
creting gas mainly by filamentary flows. The accretion pro- 
cess heats up the gas till T>7 x 10'^ K and accelerates it 
developing super sonic shocks. Near z ~ 12, at the densest 
main halo central regions the post shocked primordial gas 
has created a lot of both H2 and HD molecules in a very effi- 
cient way. Furthermore, due to the compressional turbulent 
motions triggered by the roughly radial collapse, the molec- 
ular formation is enhanced as well. Because of the molecular 
formation process the turbulent gas can reach temperatures 
even below 100 K in some regions and it is able to form po- 
tential star formation regions of dense gas (n> lO^cm"^) at 
low temperatures (T ~ few x lOOK). 

Figures [1] [2] and [3] show the z projection of the gas num- 
ber density, gas temperature and gas vorticity inside a cubic 
box of ~3 kpc of side (larger than the highest refinement 
level box of ~ 2 kpc) as a function of redshift. The horizon- 
tal and vertical lines in the vorticity map are features due 
to the lower resolution outside the highest refinement inner 
box (~2 kpc of box side). 

The analysis presented below is centered on both the 
chemo-thermal and dynamical properties of the primordial 
gas. In this paper we do not analyze either the gas fragmen- 
tation process or the primordial gas clumps physical prop- 
erties. This kind of study is postponed for a future work. 



3.1 Chemistry and Turbulence 

It is well known that below ~ lO^K the primordial gas can be 
cooled down only by molecular line emissions. In this sense 
the molecules are crucial for the formation of star formation 
regions of cold and dense gas. 

The H2 molecule is the most abundant molecule in pri- 



H H2 e" 



(2) 



mordial gas (Galli & Palla 19981. Due to its symmetry it 



has not a permanent dipolar moment and only quadrupolar 
rotational transitions are allowed ( [AbgraU et al.|1982[ ). 

The main path for H2 formation in primordial environ- 
ments is the two steps channel ( Peebles fc Dickel|1968 l 



The H2 molecule is able to drop the unperturbed primordial 
gas temperature tiU T ~ 2 - 3 x lO^K and tiU r> lOOK in 
perturbed regions with high electron fraction, e.g. in post 
shocked regions of turbulent gas with high Mach number or 
in relic HII regions. 

On the other hand, the HD molecule, also present in pri- 
mordial gas, has a permanent dipolar moment which allows 
faster rotational transitions. Furthermore, due to its higher 
reduced mass it has rotational energy levels lower than the 
H2 levels, thus reaching lower minimum gas temperatures 
due to rotational transitions ( Flower fc Roueff|1999 |. 

The main path for HD production in primordial gas is 
based on the H2 molecule (e.g. Flower (2000l): 



D^ 



H2 



HD + H+. 
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Due to the reasons mentioned above the HD molecule is 
able to cool the unperturbed primordial gas below lOOK but 
under perturbed conditions the gas could reach the CMB 
temperature fioor due to the HD cooling. 



3.1.1 H2 and HD molecules 

It is interesting to understand how the gas dynamic afl'ects 
the molecular formation process in primordial environments 
and which are the consequences of this highly non linear 
process. 

Because H2 production depends on the abundance of 
free electrons (and H^) and the HD production depends 
on the H2 (and D"*") there are two paths to increase the 
molecular production rate: i) to increase the local gas den- 
sity and enhance the number of reactions, which can be a 
natural consequence in highly turbulent environments, or ii) 
to increase the number of free electrons by a given physical 
process, e.g. a ionization front. 

The first path can be the result of either a gravitational 
collapse process or a highly turbulent environment, as men- 
tioned above. An example for the first case is the collapse of 
minihalos at high redshift ( McGreer fc Bryan|2008 l; in this 
case the density enhancement due to the gravitational col- 
lapse produces enough molecules to reduce the primordial 
gas temperature to few x 100 K. An example for the sec- 
ond case can be the interaction of gas fiows onto a (massive 
enough) DM halo; the fiows can produce a turbulent enough 
environment where over-densities fill the interaction region 
( Heitsch et al.|2008 L 

In primordial environments the second path can be the 
result of post shock waves formed in the virialization process 
of massive enough DM halos. Despite the fact that previ- 
ous works based on simplified unidimensional models have 
shown that post shock waves created through the baryonic 
collapse onto DM halos with masses M< 10* Mq cannot cre- 
ate enough free electrons in order to enhance the molecular 
production ( [Johnson fc Bromm|2006] ), this is not completely 
true due to the real non uniform baryonic velocity field, i.e. 
for a given DM average circular velocity the gas content of 
the halo is able to reach velocities higher (and lower) than 
its average and then it could create shocks able to increase 
the free electron fraction even if the average circular ve- 
locity is not high enough. Furthermore, in the hierarchical 
paradigm of structure formation it is natural that the small- 
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Figure 1. The z projection of tiie mass weiglited gas number density at z=15 (top left), 13 (top riglit), 12 (bottom left) and 10 (bottom 
right). Top left: At z = 15, near the main halo central region, appear three interacting mini halos; the nearest two over-densities will 
merge. Top right: At z = 13, after the merger a spinning over-density crosses the main halo central region as a super sonic bullet; the 
over-density leaves a post shock region where molecular coolants are formed very efficiently due to both the turbulent motions and 
the enhancement in the free electron fraction. Bottom left: At 2 = 12 the spinning over-density is leaving the main halo central region 
which starts to develop more over-densities due to the gas accretion process. Bottom right: At z = 10 the gas has developed a number 
of over-densities near the main halo central region as a consequence of the combined effect of molecular cooling and gas dynamics in 
turbulent compressed regions. 



est halos interact with the biggest ones producing a merger 
process which can create shocks waves as well and conse- 
quently increase the free electron fraction due to a violent 
merger process. 

The number density vs. gas temperature plane of our 
simulation at z—15 (top left), 13 (top right), 12 (bottom 
left) and 10 (bottom right) is shown in figure [4] The fig- 
ure shows number densities n > 10cm~^ corresponding to 
the densest gas of the simulation central regions. At z—15, 
before the merger, the gas temperature is above ~400 K 
in the densest regions associated to the minihalos. After the 
merger, the gas develops very high temperature regions with 
T>10'* K due to the violent merger process. At z—13 the 
high temperature gas has T< 10'' K because of the hydro- 
gen recombination lines cooling; This temperature is still 



enough in order to both dissociate the molecular coolants 
and enhance the ionization fraction. At this redshift, the 
high density regions, n > few x 10^ cm"'', with high tem- 
perature, r< 10*, are associated to the resultant hot-dense 
spinning over density going away the main halo central re- 
gion, while the lower density hot regions correspond to gas 
expelled out by the explosive merger. At the same time, due 
to the enhanced ionization fraction the post shocked regions 
have had enough time to create HD molecules in order to 
drop the gas temperature till the CMB limit. (In our previ- 
ous study of first galaxy formation in Prieto et al. ( 2011 1 the 



mass resolution was a factor of 8 lower than in this simula- 
tion. It is well known that in order to see the HD effect on the 
gas temperature it is needed a density threshold - |McGreerl 
& Bryan (2008I-. Due to the higher resolution of this simu- 



lation here it is possible for the HD molecule dominate the 
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Figure 2. The z projection of the mass weighted gas temperature. The different maps correspond to the same redshifts of figure[T] The 
time evolution shows how after the merger a cold tail is formed behind the spinning over-density. This remaining turbulent cold gas is 
fed by the accretion near the main halo central region and is able to develop dense and cold over-densities. 



gas cooling in some regions, including the formed gas clumps 
not analyzed in this paper.) At later times, when the dense- 
hot spinning over density is out the main halo central region 
(and certainly not captured in the bottom panels of figure 
[4|, a large amount of gas reaches temperatures ~ 100 K 
in the density range 10 < n/cm~^ < 10'': potential places for 
star formation. These final conditions are the result of a 
post-shocked turbulent evolution followed by gas accretion 
onto the central region. 

At the first glance our simulation seems contradict the 



result of Shchekinov & Vasiliev (20061. These authors found 



that in order to ionize the primordial gas in a halo merger 
processes, the halo mass should be 10^ Mq, whereas in 
our simulation the merged halos have masses ~ 10'' Mq. The 
explanation for this apparent contradiction is that the result 
of Shchekinov & Vasiliev {2006) are based on an isolated 
halo-halo interaction, whereas in our case the mini halos 
collide after to be accelerated by the fewx 10^ M© main halo, 
thus the merger is much violent than an isolated ~ 10^ Mq 
halo-halo interaction. 



Figure [S] shows the number density - molecular mass 
fraction plane for both H2 (top) and HD (bottom) molecules 
at the same four different redshift as in figure [l] As men- 
tioned above, before the merging there are a number of 
regions with temperature r>400 K and number densities 
n > 10"^ cm~^ associated with the minihalos (consistent with 
figure [4|. After the merger, both figures show that a lot of 
molecules are destroyed and the molecular mass fraction is 
reduced in a number of regions (top right panels). At the 
same time appear regions where molecules have continued 
working to reduce the gas temperature allowing the forma- 
tion of denser zones. Approximately 40 Myr later (at z=12) 
the gas has developed very dense regions with n > lO'^cm"''. 
Furthermore, there are new regions with n>10^cm~^ at 
lower molecular abundance. A hundred of million of years 
later there seems to be two populations of dense regions 
with n>few x 10''cm~^: one with temperature below 300 
K and H2 mass fraction fna ^ lO"'^ and HD mass fraction 
fnD ^ 10~^ and another one with temperature above 300 K 
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Figure 3. The z projection of the z component of gas vorticity (V X v)^. The different maps correspond to the same redshifts of figure 
^ The maps show how a turbulent environment is developed by both the merger an the accretion process near the main halo central 
region. 



and ^ 2 X 10~^ and fnD > 10"''. These two different re- 
gions can be the result of two different evolutions near the 
main halo central region. The coolest one seems to be the re- 
sult of evolution onto a post shocked environment where the 
over-densities grow up in regions with an enhanced initial 
molecular mass fraction. The other over-densities seem to 
develop in regions not too much affected by the merger pro- 
cess presenting a molecular mass fraction one order of mag- 
nitude below ( McGreer fc Bryan||2008 Johnson & Bromm 
20061). 



3.1.2 Molecular cooling and gas velocity 

The figures shown in the previous section suggest that the 
primordial gas dynamic triggers the formation of abundant 
H2 and HD molecules through both the halo merging and 
the baryonic collapse process. Consequently, these molecules 
are able to cool the gas allowing it to reach T < 300K at 
density peaks, reaching temperatures even below lOOK in 
some regions. 



Based on this dynamical dependence of molecular for- 
mation it is interesting to know how the velocity modu- 
lus/Mach number is related with the number density, gas 
temperature and molecular mass fractions. Figure |6] shows 
the molecular mass fraction - velocity modulus planes for 
both the H2 (top) and the HD (bottom) molecule at the 
same redshifts of previous figures. The velocity modulus is 

computed as the value of [(llgrid.x — Wav,x)^ + (l'grid,y— «av,y)^ + 

("Ugrid.z — Wav.z)^]^''^ for cach grid; here Uav.i is the average gas 
velocity inside the analyzed box of ~ 22 comoving kpc, i.e. 
(y^ ^Erid,i)/Ai'grid, where i — x,y,z; Vgrid,i is the local grid 
velocity and A'^grid is the total number of grids inside the 
highest resolution volume. We emphasize that the analyzed 
volume is not the entire 1 Mpc^ simulation comoving box 
but the innermost highest resolution volume of ~ 2 proper 
kpc at z=10 (~ 1.4 proper kpc at z=15). Keeping this in 
mind, the gas velocity definition of «gas = «grid — i^av has all 
sense because it represents the local grid velocity without 
take into account the gas average movement of the analyzed 
volume. (The computed rms physical velocity through the 
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Figure 4. Number density - gas temperature plane at z=15 (top left), 13 (top right), 12 (bottom left) and 10 (bottom right) for number 
densities n > 10cm~^. From the top right figure is clear how the merger process disturbs the densest main halo central region creating 
many very low temperature regions. These regions seem to be related to the post shock tail behind the spinning over-density. At later 
redshifts the temperature reaches values T < 100 K in the density range 10 < n/cm~^ < 10^ which can be associated to the HD molecular 
cooling of the turbulent main halo central region. 



evolution are 8.18 km/s at z—15, 7.06 km/s at z=13, 6.21 
km/s at z=12 and 4.18 km/s at z—10.) Despite of the main 
halo circular velocity Kirc ~ 14.6 km/s, the main halo cen- 
tral regions develop gas velocities above 22 km/s, enough to 
enhance the ionization fraction and produce shock waves[^or 
at least enough to compress significantly the primordial gas 
and accelerate the molecular formation process. At z=13, 
about 40 Myr after the merger process, the gas reaches ve- 
locities near 60 km/s evidencing the violence of the mini halo 
interaction. The interaction is able to enhance the tempera- 
ture till ~ 10" K, slightly increase the free electron fraction 
and trigger an efficient molecular formation process. 

The molecular mass fraction - Mach number planes for 
both H2 (top) and HD (bottom) molecules at the same red- 
shifts as in previous figures are shown in figure [7] The Mach 
numbers were computed based on the velocity modulus de- 
fined above. From these figures we can see that through the 
evolution almost all primordial gas with n > 10 cm~^ de- 
velops super sonic waves. After the merger there are regions 
with Mach number near 50 evidencing the violence of the 
process again. At later times the Mach numbers (and veloc- 
ities) decay at the main halo central region to increase again 
at z=10, showing how after the merger the gas relaxes and 

^ For comparison the average circular velocity of a 10**Mq DM 
halo collapsed at 2 a: 10 is Dcirc ~ 21.8 km/s. 



reaches Mach numbers below 10 at the densest regions but it 
increase again promoted by the baryonic accretion process. 
These two figures show that the coolest gas has the highest 
molecular mass fraction but interestingly this coolest gas 
occupies the high velocity/Mach number regions as a proof 
that the highly supersonic turbulent gas is composed by a 
high molecular mass fraction. 

At the moment we have a picture in which the mergers 
and the gas collapse process onto the main halo triggers ef- 
ficient molecular formation in high velocity/supersonic gas 
regions allowing the gas to reach temperatures even below 
100 K. Figure [S] shows the gas temperature - velocity mod- 
ulus plane (top) and the gas temperature - Mach number 
plane (bottom), for number density above 10 cm~'^, both 
at the same redshift as before. From this figure it is clear 
how the merger process creates high density regions with 
temperature ~ 10'' K due to its violence. At lower redshifts 
this figure shows that the high density regions with n > 10^ 
cm""^ are associated with both high Mach numbers and high 
velocities suggesting that the collapsed gas converging in 
the main halo central region is compressed at velocities 10 
km/s < Urms ^ 30 km/s, which can help to increase both the 
molecular abundances and consequently the gas cooling. In 
this way the mergers and the baryonic matter collapse pro- 
cess would favor the formation of cold - dense primordial 
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Figure 5. Top: Number density - H2 mass fraction plane at tiie same redshifts of figure |4] In brown all grid cells with temperature 
in the range lOOK to 300K, red 300K< T <500K, yellow 500K< T <1000K, green 1000K< T <5000K and blue 5000K< T <10000K. 
Due to the post shock conditions the H2 molecule is able to reach very high mass fractions, //jj 0.01, which enables the gas to drop 
temperature T<200K. Bottom: Same as top but for the HD molecule. Again, due to the post-shock conditions, the HD molecule can 
be produced very efficiently and it reaches fno < 10~^. This high mass fraction allows the gas to reach temperatures r< lOOK. 
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Figure 6. Top: H2 mass fraction - velocity modulus plane at z=15 (top left), z=13 (top right), z=12 (bottom left) and z=10 (bottom 
right). Bottom: Same as top but for HD molecule. The colors show the same range of temperature as in figure [s] Both panels show 
how the merger process increases the gas temperature and creates very high velocity gas which is able to create shocks and enhance the 
free electron fraction favoring the molecular formation in posts hock regions. After the merger, the gas of the main halo central region 
decreases its velocity but it is increased again by the baryonic accretion gas process at lower redshift. This accretion process is able to 
create high velocity gas waves which compress it and increase the molecular mass fraction. 
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Figure 7. Top: H2 mass fraction - Mach number plane at the same redshifts as previous figure. Bottom: Same as top but for HD 
molecule. The colors show the same range of temperature as in figure [s] Both planes shows that Mach numbers above 10 are reached 
in regions with both high molecular abundances and low temperatures. In general the coolest regions tend to show the highest Mach 
numbers which relates turbulent regions with high molecular abundances evidencing the effect of merger and supersonic turbulence on 
the primordial gas, i.e to increase the molecular mass fraction inducing efficient cooling. 
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Figure 8. Top: Gas temperature - velocity modulus plane for n > 10 cm ^. In yellow 10^ < n/cm ^ < 10^; in green lO'^ < n/cm ^ < 
10^; in blue 10^ < n/cm-^ < 10^ and in red IQ-* < n/cm-^ < 10^. Top: The fi gure show that high density gas with n > 100 cm is able 
to reach velocities above Vrms ~ 20 km/s, high enough to enhance the free electron number density and compress the gas increasing the 
molecular mass fraction. Bottom: Same as top but for the gas temperature - Mach number relation. High density gas with temperature 
below 10"^ K shows supersonic velocities with the highest Mach numbers associated to the lowest temperature regions. The high density 
regions with n > 100 cm~^ clearly present super-sonic turbulence. This feature relates low temperature - high Mach number with high 
density reeions showine a clear relation amone suoer sonic turbulence, molecular cooline and hieh densitv reeions. 
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gas regions and creates potential primordial star formation 
clouds. 



3.2 Probability Distribution Function (PDF) 

In the hydrodynamical context, the density PDF p gives the 
probability to find a fluid element with a density between 
p and p + dp inside the analyzed volume. In other words it 
quantifies the fraction of the gas in a given range of density 
as p(p')dp'. 



Vazquez-Semadeni ( 1994 1 shows that in the high Mach 



number limit, without viscosity and neglecting the gravita- 
tional term, the normalized Euler equations for hydrody- 
namics describe a scale invariant pressureless fiuid, i.e. the 
probability to create a relative density fiuctuation in a given 
region is independent of the local density. As explained in 
|Slyz et al.|p005l ), if we assume that in a full developed tur- 
bulent fluid the density field is a random variable and the 
event of successive density increments are independent, then 
the central limit theorem states that the density PDF is log- 
normaQ 

The argument exposed above has been extensively stud- 
ied using hydrodynamical simulations of isothermal super- 



sonic turbulence without gravity (e.g. Kritsuk et al. (20071; 



Federrath et al. (20081). These simulations have shown that 



the density PDF is well represented by a lognormal distri- 
bution: 



p(ln p)d Inp = 



\/27ro-2 



1 f Inp-lnp V 



din p, 



where 



In p = 



(4) 



(5) 



with a the standard deviation of the logarithm of the den- 
sity. For a formal proof of the lognormal density PDF dis- 



tribution based on [Pope &: Ching ( 1993 1 see Nordlund & 
|Padoan| ( p99t . 

The lognormality of the PDF is valid for isothermal tur- 
bulent fluids without gravity but under realistic conditions 
the asumption of isothermality is not well justified due to the 
cooling/heating processes triggered by both the chemistry 
and the gas dynamic. Furthermore, in less than a dynam- 
ical time the selfgravity becomes unavoidable at the high 
density regions in a realistic scenario. Therefore, we do not 
expect a priori to find a perfect log-normal density PDF in 
our non-isothermal simulation with self gravity. 

Figure [9] shows the PDF of both the gas density (top) 
and the gas temperature (bottom) for redshifts from 15 
to 10. The PDFs (for density, temperature and velocity) 
were computed as the histogram of the given quantity in- 
side equal-size bins in logarithmic scale. The density PDF 
shows a clear deviation from a log-normal distribution which 
as mentioned above is a consequence of the non isothermal 
evolution with self gravity. 

It is interesting to note how the high density tail of 

^ The PDF is lognormal because the random variable density is 
the product of independent random variables instead the sum of 
random variables, in the last case the PDF should be a normal 
distribution. 



the PDF evolves with redshift; it takes higher values at 
lower redshift showing how the gas develops high density 
(n > lO^cm"^) regions through the collapse process onto the 
main DM halo; this high density regions are less than ~ 5% 
of the volume at all redshifts. 

The features of the PDF high density tail have been re- 
cently studied by [Kritsuk et al.] ( |2011[ ) in isothermal super- 
sonic turbulence simulations with selfgravity in a periodic 
box of 5 pc size. Despite of the difference in scale among 
this two simulations, the PDFs are roughly similar in shape. 
The yellow line in flgure|9]shows ap(n) oc n~^'^ whereas they 
found a n~^'^^^ power law. This behavior is explained as the 
result of the first stage of a nearby isothermal gravitational 
collapse, as mentioned in their paper. The blue line shows a 
p{n) oc (they found ap(n) oc n~^''^^^) which is explained 
as the region in which their gas cores formed by gravitational 
collapse are rotationally stable. Interestingly, the same ex- 
planations could be valid in our simulation because through 
the gas collapse process onto the main DM halo the gas 
develops over-densities at the turbulent super sonic regions 
which eventually end as primordial gas clumps of few x 100 
K due to molecular cooling. These over-densities show in- 
trinsic rotation and develop gravitationally supported disk 
structures at lower redshifts. 

The average density inside the analyzed box for all red- 
shift of flgure[9|is n ^ (1-2) x 10"^ cm"^'. The density PDF 
shows a peak at very low densities (fewx 10~^cm~"^) instead 
that at the average density like in the lognormal isothermal 
distribution. This peak is easily explained in the collapsing 
scenario: through the collapse process, the baryonic matter 
piles up near the main halo central region. While in this 
region the gas develops over-densities in a tiny percentage 
of the volume, the rest of the volume is dominated by low 
density regions due to the converging gas flow onto the main 
halo. This dynamical behavior creates more low density re- 
gions at lower redshift as shown by the density PDF and 
favors the formation of new high density regions. 

The temperature PDF in flgure|9] shows two clear peaks 
at z=10, one at ~ few K and another one at ~ few X 10^ K. 
These two peaks are not easily recognized at 2 > 13 and they 
appear at lower redshift as a consequence of the baryonic 
matter accretion process. The two peaks feature suggests the 
existence of two preferred gas states of high and low temper- 
ature co-existing inside the analyzed volume. Again, in the 
collapsing scenario the low temperature gas can be associ- 
ated to decreasing (expanded and very low) density regions 
-which are not in equilibrium with the CMB photons- due 
to the infalling gas flow onto the main halo central region. 
These regions contain the gas inside the voids surrounded 
the high density structures (this is the reason why they do 
not appear in the flgures shown above which show gas with 
n > 10 cm~'^). On the other hand, the high temperature 
regions can be associated with heated gas infalling onto the 
DM halo through the virialization process. Between these 
two temperature peaks is located the ~ few x 10^ K gas, 
where the molecular coolants are working allowing the for- 
mation of high density regions at lower redshift. Of course 
the molecular cooling region is a tiny percentage of the total 
~ few xlO^ K range which should be dominated by ther- 
mally unstable gas. 

Because we have shown that the gas dynamic is related 
to the cooling process in primordial environments, it is in- 
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Figure 9. Top: number density probability distribution function 
for different redshifts. In black: z=15, red z=13, green z=12 and 
blue z=10. This figure shows that the density PDF has not a 
symmetric log-normal distribution and a high percentage of the 
analyzed volume corresponds to low density, n<0.1 cm~^, re- 
gions. As explained in the text, the high density tail covers two 
regimes: the first one for 0.5 < n/cm""^ < 50 shows a n~^/^ behav- 
ior (purple line) associated to a nearly isothermal collapse process; 



the second one, at densities n > 50 i 



follows a n ^ behavior 



(light blue line) associated to the disk structures of collapsed gas 
clumps. Bottom; Same as top but for the gas temperature. At 
z=10, the temperature PDF shows two clear peaks consistent 
with a two phase fluid of low temperature regions coexisting with 
high temperature regions. These two peaks are not so clear at 
z>13 and develop through the gas collapse process. 



teresting to study the mean gas velocity PDF of the system. 
Figure [TOl shows the mean gas velocity PDF for the same red- 
shifts analyzed in previous sections. The two vertical lines 
are the the main halo average circular velocity (14.6 km/s 
in purple) and the average circular velocity of a 10* Mq 
halo (21.8 km/s in light blue). At z=15, before the minihalo 
merger the gas velocity is below ~20 km/s, a velocity too low 
in order to create shock waves. At z=13, after the merger, 
the gas develops regions with velocities near 60 km/s and it 



Figure 10. Gas velocity PDF at different redshifts. In black: 
z=15, red z=13, green z=12 and blue z=10. The vertical lines 
show the main halo average circular velocity of 14.6 km/s (in 
purple) and the average circular velocity of a 10* Mq halo (in 
light blue): 21.8 km/s. The time evolution of the PDF shows how 
the minihalo merger process creates regions with velocities near 
60 km/s. After that the maximum gas velocity decreases to values 
near 30 km/s mainly supported by the gas accretion process onto 
the main halo. 



shows a longer high velocity tail. At this redshift, above 1% 
of the gas has mean velocities above 21.8 km/s. At lower red- 
shifts this percentage slightly increases but does not reach 
the 2% of the volume. Furthermore, below z=13 the max- 
imum velocity slightly surpass the 30 km/s, showing that 
the very high velocities triggered by the merger are replaced 
by moderately high velocities promoted by the gas accretion 
process on the main halo. 



3.3 Structure Functions 

In order to characterize the turbulence, we computed the 
second order velocity structure function. The velocity struc- 
ture function of order p is defined as: 



Sp{l) = (|i;(x + £)- vi^)]") oc f^'*'). 



(6) 



where the velocity component v is parallel (longitudinal 
structure function) or perpendicular (transversal structure 
function) to the vector i. The spatial average is over all val- 
ues of the position x and ^(p) is the exponent of a power 
law fit to the structure function. The second order longitu- 
dinal structure function averaged over the central region of 
the simulation from z=15 to z=10 is plotted in figure [TT| It 
is well approximated by a power law, i^''^'^, in the range of 
scales I = 50 — 200 pc (light blue line) for the given range 
of redshifts. 

Taking into account that after the merger the dynamics 
of the system is supported by the matter accretion process at 
scales of Lmax ^ 500 pc, and the molecular cooling process is 
efficient (in high density regions) at scales hmim ^ few x 10 
pc, the range in which the power law is valid corresponds 
roughly to the inertial range of the system. 




Figure 11. Root square of the second order longitudinal veloc- 
ity structure function for different redshifts. In black: z=15, red 
z=13, green z=12 and blue z=10. The light blue line shows a 
52(Qi/2 oc which is a good approximation for all curves in 

the range 50pc< I < 200pc. This range of scales roughly coincides 
with the molecular cooling (energy dissipation) scale and gas ac- 
cretion (energy injection) scale, in other words the power law is 
satisfied roughly in the inertial range of the system. The pur- 
ple line correspond to the size- velocity relation found by [Larson] 
( |1981[ l in local galactic molecular clouds. 



In order to relate this result with the observed prop- 
erties of local galactic star-forming regions, where the tur- 
bulence is driven mainly by SNe's shock waves, figure [Tl] 
actually shows the square root of the second order structure 
function next to the velocity scaling law found in galactic 
molecular clouds (|Larson||1979 19811, which follow the re- 



lation S2{£) oc i'^'"'. Interestingly, on the scale of approxi- 
mately 100 pc, the velocity dispersion in the primordial gas 
of our simulation is comparable to that of the molecular gas 
in our galaxy. However, due to the larger temperature in the 
primordial gas, the Mach number of the turbulence is a few 
times smaller than in nearby molecular clouds at the same 
scale. 



3.4 Fourier Analysis 

3.4-1 Power Spectrum (PS) 

After we have shown that supersonic turbulence, triggered 
by both merger processes and gas accretion, promotes molec- 
ular formation and gas cooling in primordial environments, 
it is interesting to study how the kinetic gas energy is dis- 
tributed in their different Fourier modes. 
The kinetic energy PS is defined as: 



(7) 



where Vk is the gas velocity Fourier transform, vl its complex 
conjugate and fc = |fe| the modulus of the wave number. Ek 
represents the kinetic energy associated to Fourier modes of 
wave number between k and k + dk. The kinetic energy spec- 
trum is useful to distinguish between subsonic Kolmogorov 



Figure 12. Velocity power spectra for different redshifts. In 
black: z=15, red z=13, green z=12 and blue z=10. The purple 
line correspond to a power spectra Ek oc . This line is a good 
approximation at z=15. At lower redshift the power spectra is a 
bit steeper and it is close to a value of 2.2. The near to Burgers 
behavior is consistent with a supersonic environment. 



turbulence, with Ek oc k~^^^, and shock dominated Burgers 
turbulence, with Ek oc k~^ . 

The energy PS for our analyzed volume is shown in fig- 
ure [12] This figure shows the kinetic energy PS from z=15 
to z=10. The purple line shows an Ek oc k~^ Burgers PS. 
It is clear that at z=15 (black line) the gas energy spec- 
trum follows a nearly Burgers behavior from fc « 10 kpc~^ 
(physical scale ~600 pc) to fc ~ 100 kpc~^ (physical scale 
~60 pc). Actually, the power law exponent is a bit steeper 
at this redshift. At lower redshifts the power law exponent 
decreases and it reaches values near -2.2 at z=10. This re- 
sult is in good agreement with the second order longitudinal 
velocity structure function shown in the previous section: 
for a Ek oc fc~" the S2{£) oc £"~^, which for n 2 gives 
approximately S2{£) oc £^ as shown in figure 11 



The close-to-Burgers power law indicates that the tur- 
bulent motions triggered by both the merger process and the 
matter accretion process onto the main halo produce super- 
sonic velocities. This result is in agreement with our results 
of subsection [3.1. 2| where we have shown that the primor- 
dial gas develops supersonic A1 > 10 motions through the 
collapse process in high density regions, near the main halo 
central region. 

It is interesting to note that at high wave numbers the 
PS does not decay as fast as in supersonic isothermal simula- 
tions without selfgravity ( Kritsuk et al.|[2007 l. This feature 
can be understood because at small scales -of order few x 10 
pc in high density regions- the gravitational potential en- 
ergy is converted into kinetic energy which is distributed 
among the longitudinal and transverse velocity modes as we 
will show in the next subsection. In this sense, in the same 
way as the large scale accretion process drives the turbulent 
and shocked environment at the main halo central region, at 
small scales the gravity is able to increase the kinetic energy 
through a turbulent energy cascade which favors the forma- 
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tion of gravitationally bounded cold dense clumps (Slyz et 



al. 2005 Federrath et al. 20111 



As mentioned in subsection |3.1.2| the supersonic turbu- 
lence and the molecular cooling seem to work togheter devel- 
oping high density regions: the large scale collapse triggers 
the formation of supersonic turbulent regions where both 
the H2 and the HD are formed efficiently cooling down the 
gas and allowing gravity to dominate at few x 10 pc scales. 
These phenomena could rise the velocity powers at small 
scales compared with the non selfgravity case as shown in fig- 
ure [12] One interesting thing to take into account for future 
idealized primordial star formation simulations, e.g. |Clark| 
|et al.| ( |2011[ ), is the PS slope found in this simulation. Our 
« — 2 slope is different to the assumed -4 slope in |Clark et| 
al. (20111. This difference could be not negligible because 



it could imply more fragmentation at small scales, certainly 
an interesting point to study. 



3.4-2 Solenoidal v/s Longitudinal Modes 

The Ifelmholtz decomposition of the velocity field allows us 
to study the kinetic energy content in compressional modes 
{vc,k with V X Vc,k = 0) and in solenoidal modes {vs,k with 
V • Vs,k = 0). In Fourier space these two components are 
defined as: 



Vc,k = {vk ■ k)k/k'^. 
Vs,k = Vk - {vk ■ k)k/k^. 



(8) 

(9) 



Figure [Ts] shows the solenoidal to total kinetic energy 



ratio 



Rk — 



\{4:iik'^Vs,k ■ K,k) 
^{Ank^Vk ■ vl) 



(10) 



from z=15 to z=10. For redshifts below 15 and for all wave 
numbers shown in the plot the solenoidal modes contain 
more than 50% of the total kinetic energy, i.e. Rk > 0.5. 
Furthermore, it is very interesting to note that above k ~ 
30kpc~^ (physical scales below ~ 200pc) R^ naturally takes 
values between 0.65 and 0.70, near a 2/3 distribution found 
in previous works ( Federrath et al.|2011 | which arises natu- 
rally in a three dimensional space where for a given compres- 
sional direction there are two rotational directions. It seems 
that at z=15 the gas has not distributed the energy follow- 
ing an equipartition rule. May be it needs the merger and 
the gas accretion onto the main halo in order to redistribute 
its energy to get the 2/3 value, in other words, the gas needs 
to pass through chaotic processes to achieve equipartition of 
energy. 

Under these result we can say that both the minihalo 
merger process and the large scale gravitational collapse pro- 
cess triggered by the main DM halo creates a supersonic tur- 
bulent cascade which naturally distribute the kinetic energy 
in a ratio Ec/Es 0.5, where Es and Ec are the kinetic 
energy in solenoidal and compressional modes, respectively. 



4 DISCUSSION 

In this paper we have shown that the mini halo merger pro- 
cess inside a ~ 3 x 10^ Mq halo is able to form shocks waves 
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Figure 13. Solenoidal to total energy ratio ij^ for different red- 
shifts. In black: z=15, red z=13, green z=12 and blue z=10. For 
a long range of k the energy ratio takes values 0.65 ^0.7. 
These values are very close to the Rj; 2/3 found in high reso- 
lution hydrodynamical simulations and is a natural consequence 
of a system with fully developed turbulence. 



and supersonic turbulence in primordial gas. Furthermore, 
we have shown that the gas accretion process onto the main 
DM halo creates more supersonic turbulence and maintains 
the turbulence formed in the merger process. It is interest- 
ing that despite of the halo mass M < 10* Mq (the limit 



suggested by Johnson & Bromm (2006 1 ) it is possible to 
find gas elements with enough velocity to create shocks (due 
to both the mergers and the gas accretion). As mentioned 
above, it is not surprising if we think that the estimation 
of [Johnson fc Bromm| ( |2006| ) is based on the average halo 
circular velocity and do not take into account the velocity 
distribution of a more realistic system. The inclusion of a 
realistic velocity distribution would let the gas to reach ve- 
locities higher (and lower) than its average circular velocity 
and creates shocked regions. This could be possible for halos 
with masses well below M < lO^M© as has been previously 



shown by Greif et al. ( 2008 1 and by this work. But here we 



have also demonstrated, for the first time, that even very 
low mass halos (M ~ W^Mq), because of the presence of 
a 1O''M0 halo, can be perturbed enough to produce tur- 
bulence and shocks inside themselves, thus triggering more 
efficient cooling than without the merging process. The pre- 
vious fact does not contradict the results of Shchekinov fcl 



Vasiliev ( 2006 1 and Vasihev & Shchekinov ( 2008 1 whom by 



one dimensional simulations shown that in order to ionize 
the primordial gas in a halo merger process, the halo mass 
should be > 10^ Mq . Despite in our case the merged halos 
have masses ~ 10® Mq , they are accelerated by a ~ 10^ Mq 
halo and they can reach enough velocity to produce the gas 
ionization. 

The supersonic turbulence formed through this process 
favors the molecular formation by two main channels: i) the 
strong shocks of velocities above 22 km/s (due to both gas 
accretion and mini halo merger) are able to enhance the 
ionization fraction which catalyzes both the H2 and the HD 
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formation increasing their abundances in post shocked re- 
gions with the consequent gas coohng, and ii) the super- 
sonic turbulence compresses the gas and increases the lo- 
cal gas density favoring molecular formation, which in turn 
enhances the lost of energy rate in the over-densities. The 
first channel is possible because the electrons recombine in 
time scales of tree ^ 10 Myr while the molecular coolants are 
formed in scales of tmoi ^ few Myr, under post shock con- 
ditions in these environments. The second channel works on 
time scales t^oi ^ 10 Myr for over densities of n ~ 10^ cm~^. 
The combination of these two process, triggered by both halo 
merger and gas accretion, facilitates the gas to drop its tem- 
perature and to form high density regions of n>10^cm~'^ 
with temperatures below 300 K, in other words: it allows 
to create cold-dense gas regions which are potential places 
for star formation. The gas temperature reaches values be- 
low 100 K at regions with 50 < n/cm~^ < 10^ where the HD 
molecule seems to be a more efficient coolant than the H2 
molecule, at higher densities the gas cooling again is domi- 
nated by H2 allowing temperatures ~ (2 — 3) x 10^ K. 

The combined effect of gas dynamic (merger and accre- 
tion) and gas chemistry (molecular formation and cooling) 
is able to produce over-densities with high velocity and high 
Mach number. The highly supersonic densest regions of the 
main halo central region develop a density PDF with inter- 
esting features. In the low density region it shows a peak 
below the average density inside the analyzed volume. This 
peak is formed because through the gas accretion process 
a large percentage of the matter is piled up in the central 
region. While the accretion flows converge near the halo cen- 
ter, the regions surrounding it lose gas continuously increas- 
ing the fraction of low density volume at lower redshift. 

At the other extreme, the high density PDF shows two 
clear power law region. The first at 0.5 < n/cm~^ < 50 fol- 
lows a PDFcx n~^'^ which is explained as the region domi- 
nated by nearby isothermal over-densities (the over-densities 
have to scale as n oc in order to reproduce the -3/2 expo- 
nent) before collapse. At higher densities, n>50 cm"'', the 
PDF tail shows how the collapse process at large scales feeds 
the central region and produces the high density structures. 
The break down of the -3/2 slope shows the point where the 
locally collapsed over-densities start to pile up the accreted 
material and develop rotationally supported structures as 
a natural consequence of angular momentum conservation. 
Very interestingly the density PDF computed from our cos- 
mological simulation has similar features to the PDF found 
in idealized star formation simulations, e.g. |Kritsuk et al.| 
(Mil). 



We emphasize that the density, temperature and veloc- 
ity PDF can be used in order to characterize the gas con- 
ditions in future idealized primordial star formation simula- 
tions inside of similar high redshift environments . 

The previous evidence of supersonic environments is 
further supported by our Fourier analysis. The kinetic en- 
ergy spectrum of the turbulent gas follows a nearly Burg- 
ers behavior oc , a bit steeper at lower redshift reach- 
ing a minimum exponent of ~-2.2. This result is also sup- 
ported by the second order longitudinal structure function 
exponent C2(^) = (^^'"^. The power spectrum slope confirms 
that merger followed by the inhomogeneous accretion pro- 
cess onto the main halo develops a supersonic turbulent gas 
with shocks through the halo virialization process. An inter- 



esting feature of this spectrum is that toward small scales it 
does not decay as fast as in isothermal simulations with no 
selfgravity. This feature could be explained by taking into 
account the selfgravity of the densest regions: At the begin- 
ning, the kinetic energy of the system comes from large scale 
accretion processes. As mentioned above, after a shocked 
is produced, the turbulent environment of the gas dynam- 
ics facilitates the gas to reach high density-low temperature 
states. Under such conditions the gas is able to locally col- 
lapse gravitationally transforming their gravitational energy 
into kinetic energy at small scales producing the feature at 
high wave numbers in the power spectrum. 

The distribution of energy in solenoidal and compressive 
modes shows that the solenoidal modes content Rk ~ 2/3 
of the total kinetic energy. This energy distribution is the 
consequence of a random turbulent motion which for each 
compressional direction, e.g. x, presents two rotational di- 
rections, y and z: the kinetic energy follows an equipartition 
rule with the same amount of energy for each direction. This 
energy distribution is not constant neither in time nor in 
scale but it has well stated values 65 < -Rfc < 70 below ~100 
pc. 

All the Fourier information presented above can be very 
useful in order to set initial conditions for idealized primor- 
dial star formation simulations in the spirit of [Clark et al.| 
(20111. For instance, these authors take an initial velocity 
power spectrum with a slope -4. This value is much lower 
than our slope ~ —2. The inclusion of a -2 slope in the 
previous study could bring a more efficient fragmentation 
at small scales with a possible more efficient low mass pri- 
mordial star formation as found in [Clark et al.| ( |2011| | and 



Greif et al. ( 2011 1 who, using 3D numerical simulations from 



cosmological initial conditions, have shown that a multiple 
system of low mass (~ 0.1 — IOM0) primordial stars can be 
formed inside DM mini halos. 



Using a one zone model Schneider et al. (20111 have 



shown that without dust there is a critical metallicity {Z ~ 
lO"*^©) above which gas can fragment into > IOMq ob- 
ject. Furthermore, they found that there is a critical dust- 
to-gas ratio above the gas can fragment into low mass 
(~ 0.01 — IM0) objects. The existence of a critical metal- 
licity above which it is possible to form low mass star is ex- 
plained by the low-mass star-formation critical metallicity 
theory of [Bromm fc Loeb (2003|) . [Frebel et alT] (|2007| tested 
this theory with data of stars formed at the galactic halo, 
globular clusters and dwarf spheroidal galaxies reaching a 
good agreement among data and theory. Nevertheless, the 
solidness of this theory has been questioned by a recently 
found extremely metal poor star which seems to violate the 
theory ( Caffau et al.|20lT I in the sense that it do not has the 
amount of both oxygen and carbon needed to form a low- 
mass star as stated in the Bromm & Loeb ( 2003 1 theory. 



This finding makes possible the potential observation of low 
mass primordial stars as mentioned in [Johnson fc Khochfar| 
(Mil). 



The study presented in this work does not take into 
account the possible radiation background from previously 
formed pop III stars. Of course, a high enough number of 
photons in the Lyman- Werner band could avoid the gas 
fragmentation due to the molecular coolants destruction 
( [Haiman et al.|[l997| |. The result of this process could be 
non fragmented hot gas piled on the main halo center of 
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mass. This configuration could be responsible for the for- 
mation of super massive black holes at high redshifts (e.g. 
i Shang et ar]|2010| )). 

Our work either take into account the possible super- 
novae explosion from pop III stars previously formed inside 
the merged mini halos. For masses of the super novae pro- 
genitor in the range 15 — 40 M0, this kind of violent phe- 
nomena certainly would be able to expel the gas inside the 
DM mini halos, and for massive enough super novae progen- 
itors, in the range of 140 — 260 Mq, the explosion could be 
able to disr upt the gas inside ha los even more massive than 



10' Mq (Whalen et al. 20081. Because our main target 



is to try to understand the main processes involved in the 
formation of supersonic turbulence in unperturbed primor- 
dial environments we have neglected these two important 
feedbacks. 

Our cosmological initial conditions were set with ag — 
0.9. This value is larger than the most recent measure of 
erg = 0.8 which imply that the DM halos in our simulation 
collapse faster than they do in a more realistic scenario. Ac- 
tually, in a more realistic scenario the halo collapse process 
would occur at a redshift 2 « 0. 89(1-1- Zo. 9) — 1, where Zo,g is 
the collapse redshift with as = 0.9. For instance, the results 
that we show at z=10 should occur at z ~ 8.8. 



5 SUMMARY AND CONCLUSIONS 

We have performed hydrodynamical simulations of primor- 
dial gas from cosmological initial conditions in order to study 
how the DM merging process triggers the formation of tur- 
bulence in primordial environments and how the developed 
supersonic environment favors the formation of molecular 
coolants. The simulation, inside a lMpc//i box size, follows 
the evolution of a 3 x 10^ M0 from z ~ 120 to z=10 with 
a proper distance resolution Ax ~ 1.95 pc at this redshift, 
covering almost 5 orders of magnitude in distance scales. 
Our main conclusions can be summarized as follows: 



trum shows that it could favor the formation of low mass 
primordial stars. 

(v) The energy spectrum does not decay at high wave 
numbers as fast as shown in isothermal simulations with- 
out self-gravity. This behavior is explained by taking into 
account that after the accretion process enhances the cool- 
ing at supersonic turbulent regions the gas is able to form 
gravitationally unstable over densities. In this way, the local 
gravitational collapse enhances the velocity power at scales 
of few X 10 pc converting the gravitational energy into ki- 
netic energy. 

(vi) The energy in solenoidal (rotational) modes is Rk ~ 
2/3 of the total energy as in previous simulations of turbu- 
lent gas. 

This paper is part of our ongoing effort ( [Prieto et al.| 



2011 1 to produce numerical simulations that capture enough 



physics to understand the process of galaxy formation at the 
highest redshifts in order to shed light on what mechanisms 
could be dominating and shaping the future evolution of 
galaxies. We have shown here that the presence of relatively 
low mass halos (~ 10^ M©) dramatically influences its en- 
vironment by creating a turbulent ISM. This in turn makes 
the smaller mass halos (M ~ IO'^Mq), the most abundant 
ones at that redshift, to also develop turbulence in their ISM 
because of the merging process. The end result is that the 
production of coolants is enhanced, so much that even the 
HD molecule becomes an important coolant in some regions, 
thus producing regions of high density-low temperature that 
could be sites of star formation. We are now investigating 
how universal the process presented in this paper is by an- 
alyzing the whole IMpc box simulation and including the 
low mass and high mass dark matter halos. We will present 
our findings in a forthcoming paper. We can at this point 
speculate that, if this mechanism turns out to be universal 
at such high redshift, it could also potentially work at lower 
redshifts, so that the initial turbulence in the ISM of a galaxy 
could be due to the merging process, which is universal in 
the standard LCDM scenario. 



(i) Both the baryonic accretion process through filaments 
and the DM minihalo merging onto the main halo is able 
to produce a supersonic {A4 > 10) turbulent and shocked 
(«rms ^ 22 km/s) environment where the H2 and the HD 
molecules are formed efficiently. 

(ii) The non equilibrium molecular formation triggered 
by the supersonic turbulent environment creates a number 
of regions with temperature r<300 K. Some regions with 
number density 50 < n < 10^ cm~^ reach temperatures be- 
low 100 K evidencing the effect of the HD molecule as the 
main coolant. At higher densities the H2 is the more efficient 
coolant. All these low temperature-high density regions are 
potential star formation places. 

(iii) The Fourier analysis of the velocity field shows that 
the kinetic energy has a nearly Burgers spectrum oc al- 
beit a bit steeper at lower redshift with a minimum exponent 
of «-2.2. This power law confirms the previous statement 
of supersonic and shocked environment. This spectrum ap- 
proximately implies a 5*2 (^) oc second order longitudinal 
structure function, steeper than the S2{£) oc £^'''^ behavior 
observed in local molecular clouds. 

(iv) The characterization of the turbulence power spec- 
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